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Critical Nature of Non- Fermi Liquid in Spin 3/2 Multipolar Kondo Model 
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A multipolar Kondo model of an impurity spin Si = 3/2 interacting with conduction electrons 
with spin s c — 3/2 is investigated using boundary conformal field theory. A two-channel Kondo 
(2CK) -like non-Fermi liquid (NFL) under the particle-hole symmetry is derived explicitly using 
a "superspin absorption" in the sector of a hidden symmetry, SO(5). We discuss the difference 
between the usual spin- 1/2 2CK NFL fixed point and the present one. In particular, we find 
that, unlike the usual 2CK model, the low temperature impurity specific heat is proportional 
to temperature. 
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An intuitive understanding of a prototype of non- 
Fermi liquid (NFL), a multichannel Kondo model, was 
proposed by Nozieres and Blandin about 25 years ago. 1 
The understanding of the multichannel Kondo model was 
obtained by various methods, such as Bethe-ansatz, 2 ' 3 
boundary conformal field theory (BCFT), 4 ~ 8 numerical 
renormalization group (NRG), 9 ' 10 and so on. In partic- 
ular, we can compare the NRG finite size spectra and 
those of the BCFT even if we have NFL spectra and the 
agreements are excellent. 7 ' 8,11 

In the realistic model of diluted f-electron systems, a 
candidate for a two-channel Kondo (2CK) model is the 
quadrupolar Kondo model 12 for U and Ce based alloys, 
in which non-Kramers doublet in f 2 configuration plays 
an important role. Several models that have a NFL prop- 
erty were proposed, in which the Tg conduction electrons 
interact with a localized f 1 or / 2 crystalline-electric-field 
state under cubic symmetry. 13-16 These models may be 
relevant for dilute alloys of Ce 3+ or U 4+ ions, such as 
Lai_ a; Ce a ;Cu2Si2, 17 UCus-^Pda; 18 among others. In gen- 
eral, NFL behaviors are observed in real materials as 
transient phenomena, which are controlled by tempera- 
ture, pressure, alloying and so on. 

In this Letter, we reexamine an impurity Kondo model 
in which a localized f-electron with Tg symmetry inter- 
acts with Tg conduction electrons under cubic symme- 
try. This model was proposed as that of Ce^Lai-zBg. 13 
The earlier NRG calculations confirmed the existence of a 
NFL fixed point that is unstable against the particle-hole 
(PH) symmetry breaking. Indeed, Ce x L&i- x TiQ is not lo- 
cated near this NFL fixed point, so that the property of 
this NFL may not be relevant to the case of Ce^Lai-^Bg. 
However, the NFL would be realized in the system with 
small quadrupolar interactions that break PH symmetry. 
Theoretically, the origin of the NFL was thought to be 
mysterious from the conventional BCFT point of view. 
The NRG energy spectrum was similar to that of 2CK 
model but the origin was unknown. 

Although it is difficult to control PH asymmetry in 
experiments, for both the pursuit of new materials and 
theoretical interest, it is worthwhile to study the detailed 



properties of this NFL by BCFT. It is noted that a BCFT 
approach is a powerful tool for answering why NFL be- 
haviors emerge in various impurity problems. The same 
situation occurs in the case of the two-impurity Kondo 
model (the NFL is unstable against PH asymmetry of 
conduction electrons), which has been extensively inves- 
tigated so far. 19 Roughly speaking, the quantum critical 
phenomena near the antifcrromagnetic critical point at 
zero temperature 20 can be seen as a kind of a lattice gen- 
eralization of the two-impurity Kondo model (however, 
it is not so simple) . The PH asymmetry in a generalized 
2CK model was also discussed. 21 In this case, in the pres- 
ence of a double tensor interaction that breaks the PH 
symmetry, the system flows to a Fermi liquid. 

We assume the following points for the local f-electron 
state: (i) cubic symmetry, (ii) f 1 Tg ground state config- 
uration, (iii) the lowest excited state is f° and/or / 2 Ti 
singlet configuration. Under these assumptions, we can 
map the Tg index to a pseudospin-3/2 representation as 
|r 8)± |) = T (|±|) + V5|t|))/V6 and |r 8i± i) = ±|±|), 
where \j z ) is a state with the total angular momentum 
J = 5/2 and its z-component j z . 

We start with the following pseudospin-3/2 Kondo 
model (see the derivation in ref. 13) in 1-dimension con- 
sidering the s-wave scattering at the impurity site: 5 

H = H + H ^ W 

m— dip, quad, oct 

Ho = ^ J2 [ dnfifaWvix), (2) 

H m = J m J2 ^(0)(xrW.(0)-X7\(3) 

where we introduce left-moving fermion annihilation op- 
erators ipfj,(x), Fermi velocity vp and spin 3/2 dipolar 
( x dip _ B ^ quadrupolar (xj? uad = q c ) and octupolar 
(x° ct = t c ) matrices of the conduction electron (similar 
definitions for X™ of the impurity). 22 J m is the coupling 
constant of each multipoles. It is noted that under the 
assumptions (i)-(iii) the interactions are isotropic in the 
pseudospin space. 
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The "conformal embedding" often found in literatures 
is SU(4H SU(2)i ® SU(2) 4 . n The SU(2)i corresponds 
to a spin current J~(x), 23 i.e., dipole. The other SU(2) 
current, SU(2)4, is an axial charge (AC) current Q(x). 24 
In this embedding, H n can be written in the following 
Sugawara form: 



ttvf 



■H 



OO 



Q-n ■ Qn ■ + 



12 



:), (4) 



where S n and Q„ are the Fourier components of 
J(x) and Q(x), respectively. We set the system size to 
21. : A : indicates the normal ordering of the operator 
A. The energy eigenvalues of the right hand side of eq. 
(4) for primary states that form conformal towers 25 are 
given by 



E(q,j) 



q{q+l) , i(i + l) 



+ 



12 

0,1/2, 
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where (q = 0, 1/2, 1) and (j = 0, 1/2, • • • , 5). The ener- 
gies of the descendant states 25 are given by E(q,j) + n, 
where n is a positive integer and indicates the PH exci- 
tations from the ground state. 6 The descendant states 
generally have different quantum numbers from those 
of the corresponding primary state in the same confor- 
mal tower. The form of eq. (4) is suitable to the case of 
J quad — Joct = 0. In this case, the impurity spin can be 
absorbed into J~ n 4 and the predicted NFL spectra are in 
complete agreement with the NRG results. 11 ' 13,15 

In the presence of quadru- and octupolar interactions, 
the above SU(2)i absorption is not applicable. As noted 
by Wu et al. the spin 3/2 fermionic system has an exact 
SO(5) symmetry under some conditions. 26 In the follow- 
ing, we transform the Hamiltonian (1) into the SO(5) lan- 
guage, i.e., the embedding is SU(4)->SO(5) 2 «) SU(2) 4 . 

First, quadrupolar matrices of s c = 3/2 have the fol- 
lowing forms: 
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where 1 is an identity matrix and the spin 3/2 oper- 
ators (i — x,y,z) are constructed on the basis of 
H4>3,ipi,i)=±, ip=i)- We have introduced five Dirac ma- 
trices r a (1 < a < 5). From eqs. (6)-(10), we can define 
ten SO (5) generators r afc as r ab = ±[T a , T b }. These gen- 
erators satisfy the SO (5) commutation relations: 

[r ab , r cd ] = -2i(s bc r ad - 5 ac r bd - s bd r ac + 5 ad r bc ). (11) 

The important point is that i?di P and H oct are written 
by the ten r ab when the condition Jdi p = J oc t is satisfied. 
Thus, under this condition, each multipolar Hamiltonian 
can be written in the SO(5)'s vector (n") and gener- 
ators Lj b {Lf') of the impurity (conduction electron), as 
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Table I. Low-energy spectrum of free Hamiltonian for nondegen- 
erate ground state. The first and the third columns are the AC 
and SO(5) indices, respectively. The second column is the eigen- 
value of the spin current in eq. (5). E(q, u) and E(q,j) are mea- 
sured from the ground state using eqs. (17) and (5), respectively. 



follows: 



#quad = 12J quad ]T<(0K, (12) 
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The free Hamiltonian Hq in (1) is expressed using the 
SO (5) generators 



/ °° 1 1 

—H = V (-: fi_„. Q n :+-:£_„•£„ :), 

ttvp ^4 8 / 



(16) 



where (C n ) ab is the n-th Fourier component of SO(5) 
current Lf(x): C n = (£ 12 , £ 13 , • • • ,Cf,Cf) and sat- 
isfies a level 2 SO(5) Kac-Moody algebra. The energy 
eigenvalues E(q, u) for the primary states are expressed 
by the Casimir operators in each sector 



E(q,u) 



q(q + 1) £, 



(17) 



where q — 0, 1/2, 1 for the AC sector (this is as same as 
in eq. (5)), and u = 1 (identity), 4 (spinor) ,5 (vector) 
for the SO (5) sector, u indicates the dimension of the 
representation in the SO (5), see Fig. 1. There are three 
primary fields in both the AC and the SO(5) sectors. The 
values of c u are ci = 0, C4 = 5/2 and C5 = 4. The non- 
interacting energy spectra for the nondegenerate ground 
state are shown in Table I together with the indices of 
j in eq. (5). We can see that both eqs. (5) and (17) can 
reproduce the spectrum of free s c = 3/2 fermions. 

Next, we consider eqs. (12), (13) and (16) together. It 
is noted that the total Hamiltonian (1) is written as 
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Fig. 1. Hi — diagrams associated with the SO(5) multiplies 
that enter the low-energy spectrum, where 2H\ = T 15 and 
2H2 = r 23 forming the Cartan subalgebra of SO (5) group. In 
(a), the corresponding components of the spinor are indicated. 
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The values at possible fixed points, the J^[ ip 
J* ct are determined (assumed) as 

TTVp 

~20l 

The case in which J qua d 7^ is discussed later. At these 
values of the couplings, the impurity SO(5) "superspin" 
can be absorbed into the conduction electron "super- 
spin" current, generating the new SO(5) "superspin" C n . 
This impurity absorption does not affect the SO (5) Kac- 
Moody algebra, except the gluing conditions in Table I, 
which is much the same as in the case of the multichannel 
Kondo model. 4 

At this stage, there is a need to introduce a suitable 
fusion rule to generate the nontrivial spectrum. Because 
the interaction is only in the SO (5) sector, any fusion 
in the AC sector is unphysical. An important point is 
that the impurity is described in the spin 3/2 represen- 
tation (i.e., dimension 4). We find that the desired fu- 
sion is a direct product of 4 representation in the SO (5) 
sector to each state in the spectra of the free Hamil- 
tonian. To execute this, the following formulae are use- 
ful (easily deducible from Fig. 1): 1 ® 4 = 4, 4 ® 4 = 
1 8 5 ® 10 and 4 5 = 4 16. 

Thus, even after absorbing the impurity "superspin" , 
we can calculate the energy spectra at this fixed point by 
using eq. (17). The obtained spectra are shown in Table 
11(a). As shown by previous NRG studies, 13, 14 the NFL 
spectra are the same as those of the 2CK NFL if the AC 
sector in the present model and the spin sector in the 
2CK model are interchanged. Indeed, the existence of a 
SO (5) symmetry in the 2CK model has been pointed out 
before in literatures. 7 
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Table II. (a) Low energy spectrum of NFL fixed point with 
^(9i u ) 5: !• I n t ne third column, the ground state energy 3/16 
is subtracted, (b) Operator contents at the NFL fixed point. We 
show the operators with A < 1 where A means the scaling di- 
mension of the corresponding operator. 



The operator contents at the NFL fixed point are ob- 
tained by applying a double fusion 4 in the SO (5) sec- 
tor (see Table 11(b)). Again, we obtain the same scal- 
ing dimensions A as those at the 2CK fixed point with 
the above-mentioned interchange. The operator with 
(q, u) = (0,5) is the SO(5) vector 0so(5)j which cor- 
responds to the quadrupolar operator in the original 
spin 3/2 representation. This operator has A = 1/2, 
so that the corresponding local quadrupolar susceptibil- 
ity diverges logarithmically at low temperatures. This is 
consistent with the result of NRG studies. 13 ' 14 The local 
charge and pair field susceptibility would diverge at low 
temperatures, because the operator (<?, u) = (1, 1), </>aci 
which is the AC vector, has A = 1/2. The dipolar suscep- 
tibility in addition to the octupolar one, is classified in 
10 representation with A = 1. This means that the dipo- 
lar and the octupolar susceptibilities \ are n °t singular, 
which is also consistent with the NRG calculation. 13, 14 

Next we discuss the stability of the NFL fixed point 
against various perturbations. 

a) In the presence of a uniaxial distortion, the term 
fiQn'j, which breaks the SO (5) symmetry, appears in the 
effective Hamiltonian. This term allows </>so(5) with A = 
1/2 to appear. Thus, the conjugate field \iq becomes a 
relevant perturbation and the NFL fixed point becomes 
unstable. 

b) In the presence of a potential scattering at the impu- 
rity site V J2 M ^(0)^(0) = VQ z (0) + const., the SU(2) 
symmetry in the AC sector is broken. In this case 



AC' 



the component q z — of </>ac with A = 1/2, can appear 
in the effective Hamiltonian. Again, the NFL fixed point 
becomes unstable. The PH symmetry can also be broken 
by the quadrupolar interaction of eq. (12). Indeed, the 
system flows into the Fermi liquid fixed point of SU(4) 
Coqblin-Schricffcr model without PH symmetry. 13 

c) The exchange anisotropy in the SO (5) sector is ir- 
relevant, because the marginal operator (q, u) = (0, 10) 
cannot appear under the time-reversal symmetry as dis- 
cussed in Ref. 8 (that is, our assumption Jd; p = J oc t, 
does not affect the present result). 27 

In early NRG studies, 13, 14 the NFL fixed point was 
considered to be similar to the 2CK fixed point. But 
the question arises, " Is the present NFL really equiva- 
lent to the 2CK ?" To verify this, we consider the ef- 
fective Hamiltonian near the NFL fixed point. If the 
present NFL and the NFL of the 2CK were equivalent, 
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Fig. 2. Ci m p/T vs T and Imxc(iij) vs £j. V is the strength of 
potential scattering at the impurity site. The parameters used 
are J,jip = Joct = 0.XD and J qua d = 0.0-D, where D is the half 
of the bandwidth of conduction electrons. 



find that the leading irrelevant operator at the NFL fixed 
point is "Fermi liquid-like" in contrast to its 2CK-like 
NFL spectra. All predictions are consistent with earlier 
works and the present NRG calculation. In particular, 
the low temperature impurity specific heat is propor- 
tional to temperature, the Wilson ratio Rw = 8/5 and 
the local charge susceptibility of conduction electrons di- 
verges at zero temperature. These results remarkably dis- 
tinguish the present NFL from the 2CK model. 
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the leading irrelevant operator would be Q_i • </>ac with 
A = 1/2 + 1 = 3/2. This operator, however, is not physi- 
cally adequate because the impurity absorption occurs in 
the SO (5) sector. The low-energy effective Hamiltonian 
should be made by using operators in the SO(5) sector 
and these should transform as a singlet. 

Because £° 6 i0go(5) (1 — a ,b,c < 5) type operators 
with A = 3/2 cannot become a singlet, the leading ir- 
relevant operator should be the energy-momentum op- 
erator J^ab '■ L c b (Q) L f (0) : with A = 2. This leads us 
to an important conclusion: impurity specific heat does 
not diverge at low temperatures unlike in the 2CK case, 
but shows a linear temperature behavior. As we men- 
tioned above, \ lS n °t singular, so the Wilson ratio 
Rw = {8x/x)/{SC/G) is calculated using the confor- 
mal charge in each of sectors cso(5) and cac as Rw — 
(cso(5) +cac)/c S o(5) - (5/2 + 3/2)/(5/2) = 8/5(C: spe- 
cific heat). 6 This result is different from that of the 2CK 
case: Rw — 8/3. 

Finally, to check the new results above, we show the 
results of the NRG calculations of the impurity specific 
heat Ci m p and the z-component of the dynamical AC 
susceptibility for localized conduction electrons Xc{<*>) at 
zero temperature in Fig. 2. Note that x c (w) is the dy- 
namical charge susceptibility of the conduction electrons 
at the impurity site. We used logarithmic discretization 
parameter A = 3 and kept 400 states at each NRG 
step. It was confirmed that Ci mp /T is constant at low 
temperatures and that lmx(0) 7^ indicates the di- 
vergence of the conduction electron's charge susceptibil- 
ity. We can see that the divergence of the charge sus- 
ceptibility is suppressed when PH symmetry is broken 
(V case). We also obtained the residual impurity 
entropy So — ks log \/2, which is the same value as that 
at the 2CK NFL fixed point. An increase in the data 
V/D = 10~ 4 of the Ci m p/T is the natural consequence 
of the entropy release from ks log \[2 to 0. 

In summary, we have investigated a multipolar Kondo 
model with Si = 3/2 and s c = 3/2 using BCFT and 
NRG. The 2CK-like NFL fixed point observed in the 
earlier NRG calculation is explicitly derived using the 
"superspin" absorption in a hidden symmetry SO (5). We 



1) P. Nozieres and A. Blandin: J. Physique 41 (1980) 193. 

2) P. B. Wiegman and A. M. Tsvelik: JETP Lett. 38 (1983) 591. 

3) N. Andrei and C. Destri: Phys. Rev. Lett. 52 (1984) 364. 

4) I. Affleck and A. W. W. Ludwig: Nucl. Phys. B360 (1991) 641. 

5) I. Affleck and A. W. W. Ludwig: Nucl. Phys. B352 (1991) 849. 

6) I. Affleck: Nucl. Phys. B336 (1990) 517. 

7) I. Affleck, A. W. W. Ludwig, H.-B. Pang and D. L. Cox: Phys. 
Rev. B 45 (1992) 7918. 

8) I. Affleck, A. W. W. Ludwig and B. A. Jones: Phys. Rev. B 
52 (1995) 9528. 

9) D. M. Cragg, P. Lloyd and P. Nozieres: J. Phys. C 13 (1980) 
803. 

10) H.-B. Pang and D. L. Cox: Phys. Rev. B 44 (1991) 9454. 

11) T.-S. Kim, L. N. Oliveira and D. L. Cox: Phys. Rev. B 55 
(1996) 12460. 

12) D. L. Cox: Phys. Rev. Lett. 59 (1987) 1240. 

13) M. Koga: Phys. Rev. B 61 (1999) 395. 

14) H. Kusunose and Y. Kuramoto: J. Phys. Soc. Jpn. 68 (1999) 
3960. 

15) M. Koga, G. Zarand and D. L. Cox: Phys. Rev. Lett. 83 (1999) 
2421. 

16) T.-S. Kim and D. L. Cox: Phys. Rev. B 54 (1996) 6494. 

17) B. Andraka: Phys. Rev. B 49 (1994) 3589. 

18) M. C. de Andrade, R. Chau, R. P. Dickey, N. R. Dilley, E. J. 
Freeman, D. A. Gajewski, M. B. Maple, R. Movshovich, A. H. 
Castro Neto, G. Castilla and B. A. Jones: Phys. Rev. Lett. 81 
(1998) 5620. 

19) B. A. Jones and C. M. Varma: Phys. Rev. Lett. 58 (1987) 843. 

20) G. R. Stewart: Rev. Mod. Phys. 73 (2001) 797. 

21) H.-B. Pang: Phys. Rev. Lett. 73 (1994) 2736. 

22) The definitions of the present multipolar operators (Q/, Tj) 
are related to those used by R. Siina et al. (O's, T's), e.g., 
Q^ 2 -r 2 = ^o^/2 and TrJ = 2T z a . See R. Siina, H. Shiba 
and P. Thalmeier: J. Phys. Soc. Jpn. 66 (1997) 1741. 

23) The spin density operator is defined as 

24) The axial charge density operator is defined as 

Q,W = E (1 (^(*W-V2), 

Q+(x) = (Q-(a:))t = i>l /2 (x)^_ 3/2 (x) - *| /2 W^ 1/2 (x). 

25) For a review, see C. Itzykson and J.-M. Drouffe, Statistical 
Field Theory (Cambridge University Press, Cambridge, Eng- 
land, 1989), Vol 2. 

26) C. Wu, J.-P. Hu and S.-C. Zhang: Phys. Rev. Lett. 91 (2003) 
186402. 

27) In the J oc t = case, the impurity spin seems to be absorbed 
as if Si = 1/2 with an even-odd interchange of the NRG it- 
eration number, resulting SU(2) NFL fixed point. 13 Although 
the precise fusion is not clear, if we admit Sj = 1/2 fusion, 
Joct might be irrelevant at the fixed point. This might suggest 
the existence of 1st order line between two NFL fixed points. 



